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ABSTRACT 


Work  has  been  completed  on  a  generalized  method  for  the 
analytical  solution  of  complex  multielement  chemical  equi¬ 
libria  under  all  conditions  of  pressure,  temperature,  and 
rqass  balance.  The  programs  are  written  for  the  IBM  7090 
and  use  a  unique  technique  of  majors  and  minors  to  guarantee 
automatic  convergence.  Variations  of  the  basic  solution 
method  are  applied  to  problems  in  chemical  synthesis,  ther¬ 
mal  stability,  and  chemical  compatibility.  The  programs 
are  also  used  to  obtain  theoretical  rocket  prop>ellant  per¬ 
formance  and  Mollier  diagrams. 
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NOMENCLATURE 


a 
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c 

p 

C 


C 
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g 

h 
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J 
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speed  of  sovind 

number  of  atoms  of  element  in  jth  chemical  constituent 
area 

specific  heat 
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total  specific  heat  of  mixtures 

molar  heat  capacity 
thrust  coefficient 

coefficient  matrix  of  the  constituent  based  on  element 
Columns  (column  matrix) 

inverted  coefficient  matrix  of  major  products  based  on  element 
rows 

free  energy 
gravitational  constant 
specific  enthalpy 
molar  enthalpy 
heat  of  formation 

arbitrary  constant  associated  with  elements 

sum  of  sensible  energy  and  heat  of  formation 
specific  impulse 

mechanical  equivalent  of  heat 
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m  weight  flow  rate 

M  molecular  weight 

n  moles 

N  weight  of  gas  or  solid 

p  pressure 

(p  )  partial  pressure  matrix  of  major  products  based  on  element 

rows  (row  matrix) 

R  gas  constant 

S  molar  entropy 

8  specific  entropy 

T  temperature 
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X  weight  fraction  of  gases 
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Y  isentropic  exponent 
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total  number  of  atoms  of  iOi  element  in  system 
Superscripts 

g  gaseous  species 

s  condensed  species 

Subscripts 
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g  gaseous  species 
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constant  pressure 
condensed  spi-cies 
throat  conditions 


constant  volume 


I.  INTRODUCTION 


The  position  of  equilibrium  in  complex  (multielement)  chemical  systems 
cati  be  computed  from  first  principles.  Given  selected  thermodynamic  and 
tlicrmochemical  data  associated  with  reactants  and  all  possible  products,  the 
distribution  between  products  and  reactants  may  be  uniquely  determined  at  all 
pressures,  temperatures,  and  stoichiometries.  Required  data  include:  heat 
capacity,  entropy,  and  heat  of  formation  of  products  and  reactants. 

The  thermodynamic  equations  that  describe  systems  which  contain  many 
elements  cannot  be  solved  in  closed  form  and  require  iterative  procedures. 
These  procedures  use  thermodynamic  properties  of  constituents  as  input  data. 
It  is  the  purpose  of  this  report  to  outline  the  capabilities  and  applications  of  a 
unique  computational  procedure  for  the  solution  of  multielement  chemical  equi¬ 
libria.  The  approach  uses  a  technique  of  majors  and  minors  to  assure  auto¬ 
matic  convJsrgence. 
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II.  OUTLINE  OF  PROGRAM 


Computer  programs  which  solve  complex  chemical  equilibria  may  be 
divided  into  two  categories.  One  type  defines  equilibria  between  constituents 
and  their  atomic  gases:  the  other,  between  the  constituents  which  are  present 
in  major  quantity  and  the  remainder  of  the  products.  These  are  known  as 
methods  of  minors  and  majors,  respectively. 

The  former  has  the  advantage  of  being  general  in  approach,  but  it  is 
logically  complex.  It  also  leads  to  numerical  difficulties  at  low  temperatures 
where  the  concentrations  of  atomic  gases  become  vanishingly  small.  The  lat¬ 
ter  has  the  advantage  of  rapid  and  certain  convergence  provided  products  are 
selected  which  are  indeed  present  in  major  quantity'.  However,  during  expan¬ 
sion  (and  cooling)  of  equilibrium  products,  the  major  products  generally  change 
and  provision  must  be  made  to  examine  and  re  select  major  products. 

The  program  system  described  herein  combines  the  flexibility  of  the 
minors  approach  with  the  convergence  advantages  of  the  majors  technique.  The 
detailed  logic  and  numerical  methods  utilized  by  the  program  will  not  be  pre¬ 
sented  here,  but  they  will  be  the  subject  of  a  future  report.  The  calculational 
procedure  is  qualitatively  the  following: 

a)  The  problem  is  solved  at  3000  K  using  the  minors  technique.  At 
3000  K,  the  concentration  of  atomic  gases  is  adequately  high  and 
convergence  is  assured. 

Since  iterative  procedures  are  utilized,  it  is  necessary  to  initiate 
the  problem  with  an  estimate  of  the  concentrations  of  atomic  gases. 
These  estimates  are  stored  internally  within  the  program  and  ob¬ 
viate  guesses  by  the  program  user.  Mass  balance  and  pressure 
balance  equations  are  applied,  and  the  solution  at  3000  K  is  attained 
even  though  the  initial  atomic  gas  concentration  guesses  may  be 
incorrect  by  factors  of  10*^. 
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I  ( 

!  b)  When  initial  converRonci!  has  hfen  altaincd,  a  ('cnc?  ra  I i/arrl  major 

product  stilection  roiitin*^  is  conirol.  A  set  of  major 

products  may  be  defined  as  lliat  list  of  products  in  which  each 
element  appears  at  least  once  and  which  contains  the  majority  of 
the  elements.^  As  many  major  products  are  chosen  as  there  are 
elemental  types.  Only  gaseous  products  are  considered  as  candi¬ 
dates.  The  routine  can  select  a  new  set  of  major  products  at  any 
time  during  the  problem  solution. 

A.  Basic  Equations 

Some  of  the  basic  equations  which  are  solved  in  the  program  may  be 
represented  by  the  following: 

1 .  Conservation  of  Atomic  Types  (Mass  Balance) 

(1) 


(2) 


(3) 


(T .  =  X  -  3..  .  n. 


2.  Conservation  of  Pressure 


j=l 


3.  Chemical  Equilibrium 


K. 


J  n 

T7‘(p.)a.  . 
1 


The  mathematically  rigorous  definition  will  be  presented  in  the  previously 
mentioned,  future  report. 
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or  using  matrix  notation  for  thf^  gcMioralized  major  product  solution 


In  =  >«  Pj  -  (C.M.)"*  (C^)  (P^) 


and,  since 


AF  =  ah  -  TAS  =  -RT  In  K 


we  have 


In  K. 
3 


-AH 

"  “RT 


+ 


AS 

"IT 


4.  Conservation  of  Static  Enthalpy 

P 

h  =5^  nj(H®)j 

3=1 


(4) 


(5) 


(6) 


B.  Minors  Technique 

Using  the  method  of  minors,  equilibrium  constants  for  constituents  are 
computed  in  terms  of  equilibrium  between  the  constituents  and  their  atomic  gases. 
To  illustrate,  consider  the  system  H,  O,  OH,  and  H2O;  equilibria  are  set  up 
in  terms  of 

OH  =  O  +  H  (7) 

HjO  =  O  +  2H  (8) 

The  thermodynamic  state  of  the  system  is  fixed  by  specifying  the  total 
amounts  of  the  elements,  (Tj,  and  the  total  pressure,  p.  The  state  of  the 
system  can  be  redefined  in  terms  of  pressure,  temperature,  and  concentra¬ 
tions  of  atomic  gases,  from  which  the  concentration  of  the  chemical  constituents 
may  be  computed. 


a 


a 
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For  solution  at  a  given  tempi' rutiire,  estimates  of  the  concentration  of 
atomic  gases  are  made  from  which  concentration  of  constituents  follow;  the 
mass  balance  and  total  pressure  restraints  arc.  of  course,  taken  into  account. 
The  equilibrium  constants  for  OH  and  H2O  are.  respectively, 


From  thf>  In  values,  the  partial  pressures  of  the  jth  chemical 
constituents  are  computed  using  Eq.  (3).  Where  the  problem  is  one  of  com¬ 
bustion  at  constant  pressure,  and  a  flame  temperature  is  sought,  the  restraint 
of  Eq.  (6)  must  be  considered. 

Estimates  of  concentration  (and  temperature,  when  necessary)  are 
repeated  until  convergence  is  obtained:  the  iteration  controls  and  sequencing 
are  complex  and  are  beyond  the  scope  of  this  report. 

C .  Majors  Technique 

Each  element  in  a  problem  requires  a  major  product,  and  all  elements 
must  be  represented  in  the  major  products  at  least  once.  If  we  select  H2O 
and  H  as  the  major  products  in  the  system  H,  O,  OH,  and  H^O,  then 

H^O  -  H  OH  (11) 
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(12) 


H2O  -  2H  =  O 


and  all  products  are  represented.  Equations  (9)  and  (10)  then  become 


•  *^OH  ^ 

+ 

^H20  "" 

+ 

For  each  nonmajor  product,  one  equilibrium  constant  needs  to  be  computed. 
D.  Program  Construction 

The  program  lu  coded  in  the  FORTRAN  system  ard  will  run  on  any  IBM 
7090  EDPM  with  two  on-line  channels  of  at  least  7-tape  units  per  channel. 

The  program  will  consider  up  to  150  chemical  constituents  and  a  maximum 
of  15  chemical  elements  per  problem. 


Each  complete  program  consists  of  an  input  or  precomputation  section, 
the  main  computation  program,  and  an  output  program.  A  schematic  is  shown 
in  Fig.  1 . 


When  the  program  recognises  the  chemical  elements  involved  in  the 
problem,  the  thermodata  tape  is  searched,  and  all  constituents  which  contain 
these  elements  are  transferred  into  core  with  their  corresponding  tabular 
thermodynamic  data.  A  program  sophistication  allows  the  program  user  to 
specify  these  constituents  which  are  to  be  considered,  if  he  so  desires. 
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COMPUTER  FILE  STORAGE 
WHEN  NOT  USED 


THERMO 

DATA 

FILE  TAPE 


THIRD  CHAIN 


PROGRAM  FILE  TAPE 


r 

IVtTEM 
SMCIFICATIOM  AND 
COMtINATION  DATA 


FIRST  CHAIN 


READS  INPUT  DATA 
CARDS.  CHECKS  FOR 
ERRORS  AND  PRODUCES 
MASS  BALANCE  DATA 


I 


ALL  PRODUCTS 


SORTS  POSSIBLE 
COMBUSTION  PRODUCTS 
FROM  MAIN  DATA  TAPE 
AND  PUTS  THEM  ON 
FIRST  BUFFER  TAPE 


X 


CALCULATES  TOTAL 
ATOMIC  MASS  AND 
ENERGY  FOR  MASS  AND 
ENERGY  CONSERVATION 
REQUIREMENTS 


T 


THERMO  BUFFER  TAPE 


STANDARD  LOOP  \ 
FOR  EACH  PROBLEM  ) 


IF  PROGRAM  ITERATION 
TEMPERATURE  MOVES 
BEYOND  INCLUDED 
CORE  RANGE- THE  NEW 
BLOCK  IS  READ  FROM 
TAPE  AND  THE  MTERNAL 
•ATA  STORAGE  UPDATED 


CONVERGES  THE  PROBLEM, 
REFERRING  TO  THE  BUFFER 
DATA  TAPE  AND  PLACING 
ANSWERS  ON  SECOND 
BUFFER  TAPE 


READS  THE  ANSWER  DATA 
BLOCKS  FROM  THE  ANSWER 
BUFFER  TAPE  AND  PROD¬ 
UCES  OUTPUT  PRINTS 


ANSWER  BUFFER 
TAPE 


FULL  ANSWER 
LISTINGS 


Fig.  1.  Program  Construction 
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The  technique  for  the  solution  of  chemical  equilibria  is  sufficiently 
general  and  independent  of  application  so  that  it  may  be  incorporated  into  a 
variety  of  special  purpose  programs.  To  date,  the  following  computation 
programs  are  included: 

1)  Chemical  Equilibria, 

2)  Adiabatic  Flame  Composition, 

3)  Mollier  Diagrams,  and 

4)  Theoretical  Rocket  Propellant  Performance. 

In  the  FORTRAN  programming  system,  all  main  programs  form  chain  links 
and  use  the  same  program  library  and  thermodynamic  data  reference  tapes. 

The  precomputation,  or  input,  program  produces  the  mass  balance  and 
data  verification  sheets.  The  input  data  consists  of  the  following  information: 

1)  tape  unit  specification  of  all  tapes  including  input,  output,  data 

( 

library,  and  program  library; 

2)  titles; 

3)  control  indicators  to  select  the  proper  computation  program; 

4)  ingredient  formulation  consisting  of  name,  elements,  and 

coefficients  contained  and  heat  of  formation  of  the  ingredient; 

5)  combination  ratio  of  ingredients:  either  moles  of  each,  mole 
per  cent,  weight  per  cent,  or  mixture  ratio; 

6)  specifically  included  or  excluded  constituents  -  optional;  and 

7)  any  specific  data  applicable  to  the  computation  chain  involved, 
such  as  enthalpy-entropy  boundaries  in  the  Mollier  diagram. 

All  control  is  internal  within  each  computational  program;  no  computer 
console  settings  are  required.  The  only  special  treatment  involved  is  the 
loading  of  the  thermodynamic  reference  tape  and  the  program  library.  All 
iteration  loops  are  checked  for  iteration  limits,  and  an  exit  procedure  is 
initiated  in  case  of  excessive  loopings. 


a 


% 
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Thi*  pronramH  named  above  usi-  a  basic  solution  block  consisting  of 
four  subroutines.  These  perform  the  functi«jns  of  examining  major  products, 
building  and  solving  the  iteration  matrices,  and  applying  the  matrix  answers. 
The  s.tme  subroutines  are  used  in  all  main  calculation  programs.  The  vary¬ 
ing  types  t)f  solutions  are  obtained  by  means  of  control  indicators  set  in  the 
main  program  and  used  in  the  basic  block.  The  basic  block  can  construct 
,»nd  s»>lve  any  <<f  the  follt>wing  t-qualion  systems; 


1)  constant  temperature  solution  of  mass  balance  and  pressure 
balance  equations; 


Z)  solution  of  mass  and  pressure  balance  equations  plus  an 
enthalpy  balance-seeks  a  flame  temperature; 

3)  system  (1)  plus  entropy  balance  (isentropic  expansion  to  assigned 
pressure); 

4)  constant  total  enthalpy  solution  and  entropy  balance  (expansion 
to  assigned  mach  number); 

“I)  constant  entropy  point  solution  (to  obtain  (9p/9p)g  for  acoustic 
velocity  determinations);  and 


6) 


constant  pressure  point  solution  (to  obtain 
determination) . 


(9n/9T)  for  C 
P  P 


The  only  program  using  all  the  available  matrices  is  the  Rocket 
Propellant  Performance  Program.  The  Chemical  Equilibria,  Mollier  Diagram, 
and  Adiabatic  Flame  Com|x»sition  Programs  use  combinations  of  (1)  and  (2). 
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III.  THERMODYNAMIC  DATA 


From  Eqs.  (9)  and  (10)  we  see  that  thermodynamic  data  required  for 
the  solution  of  equilibrium  constants  include  entropy  and  enthalpy.  Since 


C  AT 
P 


(15) 


we  find  that  is  the  sum  of  the  sensible  enthalpy  and  chemical  energy 
(heat  of  formation).  Aside  from  entropy,  we  require  heat  capacity  and  heat 
of  formation  data. 

Data  which  are  currently  available  to  the  Computer  Program  include 
the  important  products  of  the  elements  listed  in  Table  1.  While  ionic 
species  are  not  yet  included,  provision  has  been  made  for  their  incorpora¬ 
tion  (E  is  electron).  Data  on  the  thermodynamic  data  tape  are  composed  of 
heat  of  formation  and  heat  capacity,  enthalpy,  and  entropy  in  tabular  form 
from  0  to  6000  K  in  100  deg  increments.  Table  2  illustrates  a  dump  of  a 
constituent  contained  on  the  tape. 

A.  JANAF  Thermochemical  and  Thermodynamic  Data  Tables 

In  order  to  furnish  Department  of  Defense  contractors  with  consistent 
tables  of  thermodynamic  and  thermochemical  data,  the  Joint  Army-Navy- 
Air  Force  Thermochemical  Panel  was  established  in  1960.  Under 
AF  33(616)-6149  and  numerous  other  contracts,  thermodynamic  data  were 
to  be  compiled,  generated,  evaluated,  and  disseminated  in  tabular  form  and 
on  IBM  cards.  Main  contributors  include  the  Dow  Chemical  Co.  (who  issue 
the  Tables),  National  Bureau  of  Standards,  U.  S.  Bureau  of  Mines 
(Bartlesville),  and  the  University  of  Wisconsin.  Table  3  illustrates  the 
format  of  the  tables,  while  Table  4  shows  data  documentation  observed  on 
their  reverse  side. 


Table 
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Table  3 

JANAF  Thermodynamic  Data  Table 


Be 


Beryl  I ium  (Be) 


(Ideal 

Gas) 

Mol  . 

wt.  -  s 

).013 

..  -1 

/ '  " 

If.  "s 

s 

T.'K. 

c; 

§• 

-(f-MW/T 

H'-MU 

AH? 

Ar; 

*-aKr 

• 

•000 

•000 

iNOtNlTC  - 

U*«l 

TT.t*t 

TT^Ml 

INFINITI 

IM 

4«000 

*T.1M 

94*001  - 

•004 

77*790 

Ti^aia 

•  !•!•••• 

4«000 

*«•**! 

9I*000  - 

•  4M 

TI^«*T 

TI^ITT 

-  Ta.aai 

4«00i 

91*949 

97*949 

*000 

Tl.lll 

•••111 

•  90*749 

MO 
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•••IIT 

•a.iTi 
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*900 
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i*.ii* 
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1*009 
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49.009 

->  97*974 

4«000 
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99*990 

!•»•• 

Ta^wi 

40*099 

-  91.071 

TOO 
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90*709 
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1*000 

70*107 
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000 
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9*409 

Ta^aaa 
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0*040 
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4«040 
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♦{•oai 

•  T.aa* 

ITOO 

4«0M 

90*001 

90*091 

4*077 

TT^IIl 

90*144 

•  ••aai 

1*00 

4«0M 

40*110 

90*910 

9*474 

TT^eii 

94*999 

-  9*494 

1900 

4«000 

40*971 

90*901 

9*071 

70*740 

99*910 

4*094 

1000 

4«000 

40*001 

90*090 

0*400 

79*700 

90*911 

-  ♦•laT 

ITOO 

4«000 

41*109 

97*007 

0*009 

T•••l• 

!!••» 

•  9*974 

1000 

4«000 

41*477 

97*991 

7*401 

79*900 

ll•tM 

•  l•alI 

1000 

4.060 

*1.T*« 

TT^**T 

7*090 

Tmof 

99*444 

•  !•••* 

1000 

4«000 

*i««ei 

97*779 

0*499 

yu»9t 

la^ao* 

-  9*144 

1100 
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0*099 
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Table  4 

Data  Table  Documentation 


MOnLLIOH  (tm)  ( IDEAL  OAS)  NOl.  HT* 


77.241  koal. 

■ola"^ 

2M.1B  ■ 

78.25  4  0.5  keal.  uola' 

Stata  Cenflduratlon 

^S. 

»2M.16  -  ”• 

>545  eal.  dog. 

nola’^ 

Ilaetronte  Lavala  and  MiltiDlieltlaa 

1 

8 

Ji 

£^.  en**^ 

8 

1 

M 

0.0 

1 

67943 

15 

73069 

3 

21  MO 

9 

66761 

5 

73141 

1 

42S4S 

5 

69009 

3 

73520 

5 

52062 

5 

69322 

1 

73606 

1 

54677 

1 

69634 

9 

73603 

15 

56452 

5 

70606 

15 

73667 

5 

56  7n 

9 

71002 

5 

73930 

1 

SMM 

9 

71162 

3 

74071 

15 

60167 

3 

71SC0 

1 

74117 

5 

62054 

15 

71463 

9 

74163 

1 

64426 

5 

71499 

1 

74269 

15 

64507 

3 

72030 

IB 

74301 

5 

65245 

1 

72251 

5 

74416 

15 

65649 

9 

72355 

3 

74443 

5 

67226 

3 

72446 

1 

73430 

15 

TtMK  15 

73017  5 


H—t  of  yar— tlan. 

Dtrivad  tram  •  third  IM  aiMlirsla  of  tiia  vapor  proaauro  data*  soa  Ba  erpatal  for  rararofieaa> 


Waal  capacity  and  Entropy. 

Bloatronle  lavola  and  ■ultlplloltloo  froa  C.  B.  Neora.  Natl.  Bur.  Standarda  Clre.  4d7  (IMt). 


These  data  are  used  almost  exclusively.  Data  are  updated  biannually 
by  the  Panel,  and  working  meetings  are  held  triannually.  Information  on 
almost  550  solid,  liquid,  and  gaseous  constituents  are  available  over  the 
0  to  6000K range  or  where  applicable. 

While  the  tabular  data  are  currently  available  only  to  6000  K,  no 

difficulties  are  anticipated  in  extending  this  range.  The  Tables  may  simply 

be  extended  by  appropriate  calculation  or  a  simple  routine  utilised  within 

the  computational  program  which  would  extrapolate  data  to  temperatures 

greater  than  6000  K.  There  is  no  limit  to  the  maximum  temperature,  pro- 
* 

vided  ionic  species  are  considered! 

B.  Calculation  of  Thermodynamic  Data 

The  bases  for  thermodynamic  data  calculation  are  standard  to 
statistical  mechanics. 

Assignments  for  ideal  monoatomic  gases  include  translation  and 
electronic  contributions. 

For  ideal  diatomic  gases,  contributions  from  translation,  rotation, 
vibration,  and  electronic  transitions  are  included.  When  available,  spectro¬ 
scopic  constants  are  utilized;  otherwise,  appropriate  molecular  models  are 
used.  Anharmonicity  corrections  are  made  when  adequate  data  are  available. 

For  linear  polyatomic  molecules,  anharmonic  corrections  are 
neglected. 

The  rigid  rotator,  harmonic  oscillator  approximation  is  utilised  for 
nonlinear  polyatomic  molecules,  and  anharmonic  corrections  are  neglected. 

Thermodynamic  data  of  condensed  phases  rest  upon  either  measured 
or  estimated  heat  capacity  data.  Heats  of  transition  are  either  experi¬ 
mentally  determined  or  theoretically  estimated  and  are  included  in  enthalpy 
and  entropy  totals. 


The  present  computation  system  automatically  adds  a  charge  balance 
equation  when  ionic  species  are  included  in  the  problem. 


a 


% 
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c. 


Determination  of  Thermochemical  Data 


Generally,  the  most  important  thermodynamic  or  thermochemical 
datum  is  the  heat  of  formation  of  the  constituent.  Heats  of  formation  are 
directly  obtained  from  calorimetric  measurements  or  indirectly  from  bond 
energy  contributions  as  estimated  from  spectroscopic  data  or  calorimetric 
data.  Second  law  heats  are  utilized  as  indirect  estimates  when  necessary. 

D.  Thermodata  Tape  Monitor  Programs 

Programs  which  maintain  and  update  the  thermodata  tape  are  inde¬ 
pendent  of  the  calculational  program.  The  monitor  routines  will  alter  any  of 
the  data  as  given  in  Table  2  as  well  as  add  or  delete  elements  or  constituents. 
In  every  case  of  data  alteration,  the  tape  identification  ntmiber  is  updated  as 
the  first  action.  Where  different  individuals  require  their  specific  data  tapes, 
a  number  of  tapes  may  be  written. 


IV.  PROGRAM  APPLICATIONS 


The  development  of  the  capability  to  predict  chemical  equilibria  in 
systems  under  analytical  or  experimental  investigation  is  a  valuable  adjunct 
to  problems  involving  chemical  synthesis,  materials  development  and  test¬ 
ing,  laboratory  flames,  chemical  propulsion,  and  aerodynamics.  Whenever 
possible,  thermodynamic  calculationp.  should  be  applied  to  experimental  pro¬ 
grams.  Consideration  of  the  thermodynamics  of  the  situation  can  save  a 
great  deal  of  experimentation  in  terms  of  a  priori  analysis. 

A.  Chemical  Synthesis 

Experimental  attempts  to  create  new  molecules  through  synthesis  or 
high  temperature  equilibration  can  be  placed  on  a  sound  basis  by  computing 
the  equilibrium  composition  of  products  under  specified  reaction  conditions. 
For  illustrative  purposes,  consider  the  synthesis  of  nitric  oxide  (NO)  by  the 
high  temperature  equilibration  (over  a  catalyst)  of  an  arbitrarily  chosen 
equimolar  mixture  of  nitrogen  and  oxygen.  (See  Tables  5  -10.)  Input  data 
to  the  program  include  the  following  which  must  be  specified; 

1)  molar  or  weight  concentrations  of  reactants; 

2)  total  pressure(s)  at  equilibration;  and 

3)  temperature(s)  at  equilibration. 

For  an  example  (only  a  few  points  are  illustrated),  69  temperatures  from 
6000  to  300  K  were  computed  at  a  total  pressure  of  1  atm  in  2.  6  min.  Be¬ 
tween  6000  and  1000  K,  temperature  increments  were  100  K;  1000  to  500  K, 
increments  were  50  K;  and  500  to  298.  16  K,  increments  were  25  K. 

Table  5  shows  the  problem  definition  and  mass  balance  sheets  followed 
by  an  automatic  statement  of  the  products  selected  for  consideration  and  their 
heats  of  formation  (Table  6).  Tables  7-10  show  the  computed  equilibria  at 
various  temperatures.  Included  are  their  reactant  and  product  concentrations 
and  their  individual  heat  capacities  (CP),  enthalpies,  entropies,  and  free 


u 


energies  (H-TS).  System  enthalpies,  entropies,  and  heat  eupaeitiirs  (frozen 
equilibria)  are  also  tabulated,  p'rum  the  example,  it  was  conc  luded  that  a 
maximum  concentration  of  NO  lor  the  givi;n  stoichiometry  c-nd  ]>ressure 
would  be  attained  at  approxintiatcly  JSOO  K. 

It  should  be  noted  that  the  program  can  yield  adiabatic  flame  tempera¬ 
tures.  Since  the  enthalpy'"'  of  NO  at  29S.  16  K  is 

AH,  +  H° 

1.  _  t  o 


+  21.  60  +  ^(1)  (1.89)  +  (1)  (2.  06) 
30 


=  0.85  (16) 

0.85  kcal/gm,  we  would  calculate  that  the  adiabatic  decomposition  flame 
temperature  of  NO  (initially  at  298.  16  K)  is  between  2600  and  2700  K  (see 
Tables  9  and  10).  Plots  of  system  H  vs  T  over  small  temperature  intervals 
are  linear,  and  linear  interpolation  is  justified. 

B.  High  Temperature  Stability 

In  many  cases  it  is  desired  to  predict  the  thermal  stability  or  compati¬ 
bility  of  refractories  or  other  materials  at  elevated  temperatures  in  reactive  and 
inert  atmospheres.  Tables  11-14  illustrate  thermodynamic  computations  of 
the  refractory  beryllium  oxide  in  gaseous  nitrogen  and  water  vapor.  The  cal¬ 
culations  refer  to  the  stoichiometries 

Stability  -  BeO  +  N2 

Compatibility  -  BeO  +  H^O 

in  closed  systems  and  predict  the  distribution  of  products  which  would  be 

AH-  from  Table  6  and  H°  from  Table  1. 

f  o 


Bf  CONO  0.  7.5000  152.0641  10.6747  96.0400 

OEO  CONO  0.49201340  16.0000  36.7120  33.9904  >65.2593 
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Equilibration  of  BeO  and  H-O  at  3000  K 
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Table  14 

Equilibration  of  BeO  and  H^O  at  2300  K 
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observed  after  an  infinite  length  of  time.  The  tables  show  that  in  nitrogen 
condensed  BeO  decomposition  proceeds  through  the  formation  of  the  gaseous 
BeO  polymers.  In  water  vapor,  the  gaseous  hydroxides  play  a  major  role. 

The  partial  pressures  of  the  constituents  are  simply  their  mole  fractions 
multiplied  by  the  total  pressure. 

C.  Composition  of  Flames 

The  computer  programs  can  be  of  great  value  to  problems  which  in¬ 
vestigate  flame  phenomena.  Input  data  to  the  program  would  consist  of 
oxidiser  and  fuel  stoichiometry  and  ambient  pressure.  These  conditions 
would  define  a  premixed  flame  or  a  rocket  exhaust  jet.  The  actual  flame 
temperature  at  a  given  stoichiometry  and  pressure  would  easily  be  calcu¬ 
lable  from  the  thermodynamic  and  thermochemic.al  data  of  the  oxidizer  and 
fuel. 

Tables  15  and  16  list  the  flame  compositions  of  the  system 
tetrafluoroethylene -oxygen  (C2F^  +  2O2)  at  ambient  pressure.  The  calcu¬ 
lations  were  ‘made  in  support  of  a  laboratory  program  which  will  investigate 
the  flame  in  order  to  give  insight  into  the  phenomonology  of  Teflon  ablation 
in  air.  One  of  the  surprising  results  is  the  extremely  high  concentration 
of  fluorine  atoms.  For  a  given  stoichiometry  and  pressure,  this  program 
does  not  automatically  compute  the  adiabatic  flame  temperature;  these  would 
be  hand  calculated  as  shown  in  Section  IV.  A.  ,  Chemical  Synthesis.  Oxidiaer 
and  fuel  enthalpy  is  l.OO  kcal/gm  which  yields  an  adiabatic  flame  temperature 
between  2000  -  2100  K. 

Calculations,  suchasthose  illustrated  in  Tables  15  and  16,  would  appear 
to  be  an  aid  toward  the  prediction  of  the  radiative  signatures  of  rocket  ex¬ 
haust  jets.  Since  the  flame  composition  and  temperature  is  predicted,  the 
concentration  of  the  emitters  from  first  principles  has  been  established. 

These  calculations  are,  of  course,  based  on  thermodynamic  equilibrium  and 
do  not  recognise  the  existence  of  excited  emitters  or  other  nonequilibrium 
phenomena. 
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For  the  premixed  flame,  the  specification  of  oxidizer  and  fuel  stoi¬ 
chiometry  is  reasonable  since  one  expects  and  observes  predicted  flame 
temperatures,  radiation,  and  compositions.  However,  many  processes  of 
interest  to  the  aerospace  sciences,  such  as  ablation,  must  be  approximated 
by  a  diffusion  flame.  There  is  transport  of  material  through  the  boundary 
layer,  and  the  flame  stoichiometry  is  not  well  defined.  However,  flame 
temperatures  near  the  theoretical  maximum  are  often  observed  in  diffusion 
flames  although  there  is  no  true  final  flame  temperature. 

D.  Problems  Which  Contain  Many  Elements 

Tables  17  -  19  illustrate  a  problem  which  contains  six  elements.  As 
seen  from  the  Tables,  condensed  phases  may  appear  and  disappear.  Com¬ 
putation  time  is  somewhat  increased  over  that  with  fewer  elements;  there 
are  no  increased  computational  difficulties. 

E.  Mollier  Diagrams 

Isobaric  and  isothermal  curves  on  an  enthalpy -entropy  plane  for  a 
given  stoichiometry  result  in  a  Mollier  diagram.  Such  diagrams  are  of 
general  utility  in  power-plant  design  and  thermodynamic  analysis.  Mollier 
diagrams  are  usually  presented  for  single  constituents,  but  the  computer 
routines  described  in  this  report  are  general  in  nature.  Thus,  data  for  the 
Mollier  diagrams  can  be  easily  and  automatically  computed  for  systems 
which  contain  numerous  constituents. 

Printouts  from  the  current  Mollier  Diagram  Program  are  illustrated 
in  Tables  20  -  22  and  are  produced  at  various  pressure  and  temperature  inter¬ 
vals  specified  by  the  program  user.  In  practice,  points  are  generated  by 
computing  points  at  equally  spaced  temperatures  along  an  isobar.  When  com¬ 
putations  have  been  completed  for  a  temperature  range  and  pressure,  the 
program  proceeds  to  the  next  pressure.  The  following  data  are  printed  out 
at  each  pressure  and  temperature; 
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Table  18  (continued) 

Equilibration  of  C10,F  and  L>iBH.  at  IZOO  K 
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Table  22  (continued) 
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Column 


Datum 


1  Pressure 

2  Temperature 

3  Enthalpy 

4  Entropy 

5  Indicator 

6  Point  Number 

When  condensed  phases  appear,  the  program  notes  their  presence  and  iden¬ 
tifies  them  (see  Table  22). 

The  subscript  COND  refers  to  a  constituent  whose  thermodynamic  data 
over  the  range  0  to  6000  K  is  associated  with  both  solid  and  liquid  phases. 

The  data  list  the  melting  point  to  the  nearest  100  deg,  and  the  entropy,  heat 
capacity,  and  enthalpy  are  discontinuous.  System  composition  data  are  deleted 
in  order  to  reduce  the  amount  of  output  but  may  be  produced  as  an  option.  The 
temperature  interval  selected  for  each  isobar  is  identical,  facilitating  the  con¬ 
struction  of  isotherms  on  the  diagram. 

Applications  of  Mollier  diagrams  of  complex  systems  are  numerous.  In 
the  aerospace  sciences,  obvious  application  involves  determination  of  theo¬ 
retical  rocket  propellant  performance  and  the  solution  of  a  host  of  aerodynamic 
problems.  Most  important  in  the  latter  case  is  ablation  and  the  calculation  of 
the  thermodynamic  and  transport  properties  of  flowing,  reacting  gas  mixtures. 
Since  the  viscosity  and  thermal  conductivity  of  a  gas  mixture  is  of  importance, 
gas  composition  must  be  known  before  its  transport  properties  can  be  comp- 
puted.  The  ability  to  treat  flow  with  chemical  reaction  is  mandatory  for  the 
solution  of  contemporary  aerodynamic  problems. 

While  each  point  computed  by  the  Mollier  Diagram  Program  is  an 
equilibrium  one,  shifting  equilibrium  parameters  such  as  heat  capacity,  isen- 
tropic  exponent,  and  speed  of  sound  necessary  for  some  aerodynamic  calcu¬ 
lations  are  not  computed.  The  program  could  be  modified  to  compute  these 
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parameters,  but  computer  time  woulfl  be  increased  sharply.  However,  the 
following  section  will  illustrate  how  the  parameters  may  be  derived  from  the 
current  program  which  will  be  modified  slightly. 
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F.  Calculation  of  the  Velocity  of  Sound 
The  velocity  of  sound  is 
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(17) 


and  y  the  isentropic  exponent  is 
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Using  the  Bridgman  formulas,  Eq.  (18)  is 
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Since 
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differentiation  yields 


p  dv  +  V  dp 
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and  substitution  in  Eq.  (19) 
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where  c  , 
P 


the  shifting  equilibrium  heat  capacity,  is 
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For  a  nonreacting  gas  mixture  (frozen  equilibrium) 
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=  VFTJrT-;p  =  ° 

and  Eqs.  (22)  and  (17)  are  reduced  to 
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and 
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(24) 


(25) 


(26) 


Mol«Tular  weight  changes  are  primarily  caused  by  reassociation  of  high 
temperature  species;  the  degree  of  dissociation  is  related  exponentially  to 
temperature  and  is  much  less  sensitive  to  pressure.  At  static  temperatures 
l>elow.  e.g.  ,  ISOO  K,  molecular  weight  changes  become  negligible  and  Eq.  (25) 
may  be  utilized.  At  temperatures  below  approximately  1500  K,  frozen  and 
shifting  equilibrium  parameters  arc  essentially  identical.  Table  23  illustrates 
values  for  frozen  and  shifting  y  in  the  oxygen /hyd  rogen  system. 
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In  order  to  evaluate  the*  HhiflinK  equilibrium  isentrupic  exponent  and 
the  speed  of  so\md  for  areas  on  a  Mollier  diagram,  we  need  information  on 

1)  enthalpy 

2)  molecular  weight 

3)  pressure 

4)  temperature 

5)  volume 

These  will  be  tabulated  for  each  point;  for  a  given  isobar,  the  molecular 
weight,  enthalpy,  and  specific  volume  will  be  found  on  the  printout.  The 
speed  of  sound  and  isentropic  exponent  may  then  be  evaluated  graphically. 

A  section  of  a  Mollier  diagram  is  illustrated  in  Fig.  2,  and  we  wish 

to  calculate  the  isentropic  exponent  in  the  area  of  pj  and  T^  from  Eq.  (22). 

Evaluation  of  (B  lnM)/(dln  p),.  and  (8  lnM)/(dln  T)  would  be  made  over 

^  P 

the  interval  p^,  Tj  -  Tj  and  p^,  T^  -  Pj.  T^,  respectively,  as 

(A  lnM)/(Aln  T)  .  c  is  simply  read  from  the  printout  as  the  difference  in 
P  P 

h  between  Pj,  T2  -  pj.  T^.  For  the  speed  of  sound  we  need,  in  addition, 
pv  at  P|,  Tj.  Since  the  data  are  printed  out  along  an  isobar 
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Use  of  Correct  M  When  System  Contains  Condensed  Phases 
Let 


and  for  a  system  which  contains  a  condensed  phase 
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Since 
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then  Eq.  (29)  is 


V 


N 

-IP- 


N  c*  -  N 
g  P 


N  c 


g  TJ 


N  c" 
s  p 


(30) 


(31) 


(32) 


whe  re 


c  =  total  specific  heat  of  mixtures 
P 

N  =  weight  of  gas  or  solid 


X  -  weight  fraction  of  gases 
g 


M 


molecular  weight  of  mixture 


Thus,  the  M  used  in  a  two-phase  system  is  actually  M/X^. 
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Theoretical  Rocket  Propellant  Performance  Calculations 


The  solution  of  complex  chemical  ecjuilibria  is  a  fundainiMital  require- 
m«*nt  for  the  computation  of  theoretical  rockc't  propellant  performance.  Every 
parameter  of  performance  may  be  computed  from  a  thermodynamic  analysis  of 
c«»mbustion  products  at  selected  pressure  planes  within  the  thrust  chamber. 


I. 


propellant 


Performance  Parameters 

Generally,  performance  parameters  which  characte ri/a-  a  given 
system  are 

a)  specific  impulse 

b)  I  cha racte ri .st u"  \  eloi  ity 


S2 


thrust  coefficient 


d)  c  area  ratio 

Items  (a)  and  (c)  are  a  function  of  the  area  ratio,  and  it  is  usually  desired 
that  they  be  tabulated  at  numerous  pressure  ratios  in  order  to  predict  alti¬ 
tude  and  vacuum  performance.  Also,  it  is  advantageous  to  include  tabular 
material  for  the  calculation  of  performance  resulting  from  nonoptimum 
expansion. 

2.  Specific  Impulse 

Specific  impulse  is  defined  and  calculated  from  the  relation 
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Optimum  specific  impulse  is  established  when  the  static  pressure  at  the 
nozzle  exit  is  equal  to  the  ambient  pressure,  and  Eq.  (33)  is 


(34) 


When  the  thrust  chamber  is  at  a  high  altitude  and  the  ambient  pressure 
approaches  a  vacuum,  Eq.  (33)  is 
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Both  I  .  and  I  are  a  function  of  the  nozzle  expansion  ratio  (t);  the 
opt  vac 

latter  is  defined  by 
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(36) 
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Because  the  rockft  engine  in  a  boat  engine  in  which  the  kinetic  energy  of  the 
jet  is  created  at  the  expense  of  combustion  product  enthalpy,  it  can  easily  be 
shown  that 


(37) 


The  calculation  of  from  Eq.  (37)  involves  the  difference  in 

combustion  product  enthalpy  between  the  chamber  and  the  exit  pressure  plane. 
The  calculation  proceeds  on  the  basis  of  selecting  a  chamber  pressure  and 
assigning  the  ambient  pressure  to  the  pressure  at  the  exit  plane  of  the  nozzle. 
Calculations  are  made  with  the  following  assumptions: 


a)  Propellants  are  injected  into  a  combustion  chamber  at  a 
prescribed  chamber  pressure. 

b)  Propellants  are  transformed  adiabatically  and  isenthal- 
pically  into  their  combustion  products  at  the  prescribed 
chamber  pressure. 

c)  Combustion  is  complete,  and  a  mixture  results  which  can  be 
described  by  the  perfect  gas  law. 

d)  Combustion  products  enter  a  de  Laval  nozzle;  the  contraction 
ratio  (A^/Aj)  is  infinite,  and  the  velocity  of  gases  entering 
the  convergent  portion  of  the  throat  is  negligible. 

e)  Jet  velocity  at  the  throat  is  equal  to  the  local  speed  of 
sound . 

f)  Expansion  through  the  nozzle  is  adiabatic  isentropic,  friction¬ 
less,  and  one-dimensional. 

g)  Condensed  phases  exhibit  a  negligible  volume  and  are  in 
thermal  equilibrium  with  gas  particles  while  travelling  with 
their  identical  velocities. 


The  expansion  process  may  be  treated  exactly  by  the  methods  of 
classical  thermodynamics  with  the  assumptions  of  either  frozen  or  shifting 
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equilibrium.  For  the  former,  we  assume  the  composition  of  the  chamber 
combustion  products  to  be  invariant  with  position  in  the  nozsle.  With  the 
assumption  of  shifting  equilibrium,  combustion  product  composition  continu¬ 
ally  changes  and  is  required  to  reflect  the  temperature  change  during  the 
nozzle  expansion  process.  All  performance  parameters  may  be  computed  for 
the  conditions  of  frozen  and  shifting  equilibrium. 

3.  Area  Ratio  (c) 

Thrust  chamber  geometry  is  best  characterized  by  the  area 
ratio,  (  .  Utilizing  the  equation  for  the  continuity  of  mass 


puA  =  m 


(38) 


Eq.  (36)  becomes 


The  perfect  gas  law  is 


and  c  is  then 


c 


(pu)^ 


(39) 


(40) 


(41) 


€  is  seen  to  involve  the  evaluation  of  the  properties  of  combustion  products 
at  the  throat  and  exit  pressure  plane. 

Procedures  for  evaluating  throat  properties  are  somewhat  compli¬ 
cated.  since  the  throat  position  is  not  uniquely  known  but  encompasses  only 
basic  thermodynamics.  The  significance  of  Eq.  (41)  is  that  from  a  thermo¬ 
dynamic  analysis  of  combustion  products  at  selected  stations  in  the  thrust 
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chamber,  data  on  the  design  of  thrust  chambers  may  be  deduced.  Given  a 
propellant  stoichiometry  and  enthalpy,  chamber  pressure,  and  exit  pressure, 
an  optimum  impulse  can  be  computed  and  the  area  ratio  of  the  thrust  chamber 
which  is  required  to  produce  that  exit  pressure  may  also  be  computed. 

Substituting  Eqs.  (34),  (36),  and  (38)  into  Eq.  (33) 


I 

( 


(42) 


from  which  can  be  computed  the  specific  impulse  at  any  area  ratio  (via,  exit 
pressure)  and  for  all  conditions  of  nonoptimum  expansion,  p^  ^  0.  The 
special  case.  namely,  p^  =  0,  is  also  included.  The  term 

A^/m  =  €/(pu)^  can  be  expressed  in  more  tractable  terms.  From  Eqs.  (34), 
(38).  and  (40). 


Thus,  nonoptimum  specific  impulse  can  be  computed  but  without  specification 
of  €  .  However,  the  specification  of  the  area  ratio  for  optimum  or  nonoptimum 
specific  impulse  will  always  involve  the  determination  of  throat  properties 
since  t  can  only  be  derived  from  Eq.  (41). 

4.  Derived  Performance  Parameters  -  Thrust  Coefficient  and 
Characteristic  Velocity 

For  reasons  peculiar  to  the  experimental  evaluation  of  a  propel¬ 
lant  system  in  a  thrust  chamber,  it  is  convenient  to  define  parameters  which 
yield  insight  into  the  efficiencies  of  the  combustion  and  expansion  process. 
Three  parameters  associated  with  the  thrust  chamber  which  may  be  experi¬ 
mentally  measured  are  the  chamber  pressure,  throat  area,  and  propellant 
weight  flow.  As  an  excellent  approximation  for  a  given  throat  area,  the 
chamber  pressure  is  directly  proportional  to  the  propellant  weight  flow. 

Pc  a  ^  (44) 
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Note  that  the  term  p  A./rh  has  the  dimeniion  of  seconds,  which  are  those 

c  i 

of  specific  impulse;  and  we  intuitively  conclude  that  the  term  could  be  pro¬ 
portional  to  specific  impulse.  Utilising  the  acceleration  of  gravity  we  find 
the  familiar  proportionately  constant  c'*‘ 


c 


(45) 


to  have  the  dimensions  of  velocity. 

The  characteristic  velocity,  c*^,  has  appeared  quite  often  in  the 
rocket  literature  and  has  been  the  source  of  much  misunderstanding.  Much 
of  the  confusion  associated  with  the  interpretation  of  c*^  could  be  avoided  if 
its  lack  of  fundamental  significance  could  be  kept  in  mind.  Unfortunately, 
c*'  is  neither  the  velocity  of  the  jet  at  the  throat  nor  at  the  exit  plane,  and 
the  product  of  chamber  pressure  and  throat  area  has  no  particular  signifi¬ 
cance.  We  find  Eq.  (45)  to  be  a  general  equation  which  states  that  for  a 
given  throa^  area  and  chamber  pressure,  propellant  systems  with  higher  c^ 
require  less  flow  rate  to  produce  a  given  chamber  pressure.  Equations 
(44)  and  (45)  are  statements  of  a  relationship  between  chamber  pressure 
and  propellant  flow  rate  for  a  given  propellant  system.  The  value  of  the 
parameter  c''‘  is  that  while  Eq.  (45)  may  be  evaluated  experimentally,  it 
may  also  be  evaluated  theoretically  from  a  thermodynamic  analysis  of  the 
propellants,  and  the  results  can  be  compared.  Thus,  one  obtains  an  experi¬ 
mental  indication  of  whether  the  particular  design  of  hardware  (injector  or 
thrust  chamber)  permits  the  realization  of  near  theoretical  results.  Use  of 
Eqs.  (38),  (40),  and  (45)  yields 


(46) 


c*  is  a  single  valued  function  of  propellant  stoichiometry  inlet  conditions 
and  chamber  pressure  and  is  independent  of  processes  which  take  place  be¬ 
yond  the  throat.  is  a  function  of  propellant  stoichiometry  and  chamber 
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pressure  but  includes  the  contribution  of  the  expansion  process.  We  define 
the  coefficient  of  thrust 


(47) 


and 
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(48) 


We  define  c’’*  and  Cp  to  approximately  separate  the  combustion  and  expan¬ 
sion  efficiencies.  Like  and  c*^,  Cp  may  be  computed  for  both  froeen 
and  shifting  equilibrium. 

In  general,  c^,  which  is  determined  experimentally,  is  less  than 
c*  theoretical,  while  the  reverse  is  true  for  the  thrust  coefficient.  The 
reason  is  that  combustion  is  usually  not  complete  in  the  chamber  and  will 
occur  to  some  extent  in  the  nozzle.  The  result  is  a  somewhat  lower  than 
theoretical  chamber  pressure.  Hence,  the  experimental  c*^  is  lower  than 
theoretical,  while  the  enthalpy  converted  to  velocity  in  the  nozzle  is  aug¬ 
mented  by  an  additional  enthalpy- of  combustion.  It  should  be  realized  that 
the  thrust  coefficient  is  determined  indirectly;  thrust,  chamber  pressure, 
and  propellant  weight  flow  are  measured  directly,  and  the  thrust  coefficient 
is  computed. 

5.  Theoretical  Rocket  Performance  Printout 

Tables  24  -  28  illustrate  the  output  format  of  the  Shifting  Equilib¬ 
rium  Theoretical  Rocket  Performance  Programs.  Input  data,  such  as  pro¬ 
pellant  chemical  formulas  and  heats  of  formation,  are  reiterated.  Mass 
balance  sheets  are  then  produced,  and  propellant  components  are  listed  in 
weight  per  cent,  moles,  and  mole  per  cent.  A  list  of  the  combustion  products 
selected  and  their  heats  of  formation  (kcal/mole)  follows. 
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Table  26 


Table  26  (continued) 


Table  27 
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The  problem  solved  in  the  tables  illustrates  the  typical  three- 
table  shifting  equilibrium  output  format.  While  the  selection  of  the  chamber 
pressure  is  infinitely  variable,  computations  are  executed  at  24  pressure 
ratios  (Table  26).  The  initial  pressure  is  that  within  the  combustion  chamber 
whereas  the  fifth  pressure  locates  the  exact  position  and  conditions  at  the 
throat.  At  every  pressure,  the  temperature,  enthalpy,  entropy,  heat  capacity, 
isentropic  exponent  (GAMMA),  product  molecular  weights,  and  the  concentra¬ 
tion  of  combustion  products  are  printed. 

Table  27  lists  the  rocket  performance  data.  Characteristic 
velocity  is  tabulated  as  is  optimum  specific  impulse,  expansion  ratio 
(EPSILON),  and  thrust  coefficient  (CF)  as  a  function  of  pressure  ratio. 

Table  28  lists  altitude  performance  parameters  as  a  function  of 
50  unit  area  ratios  and  is  utilized  to  read  the  thrust  coefficient  (CF  EPS)  and 
specific  impulse  (I  EPS)  for  an  ambient  pressure  of  zero  as  a  functipn  of 
area  ratio.  Optimum  specific  impulse  (I  OPT)  is  tabulated  at  the  unit  area 
ratio.  The  table  affords  the  rapid  and  accurate  calculation  of  specific  im¬ 
pulse  under  any  set  of  optimum  and  nonoptimum  area  ratios  and  altitude 
conditions.  The  exact  equation  for  these  effects  is 

I  EPS  =  I  OPT  +  A  -  Bp 


p  c  C  p  c  *€ 

r  I  OPT  +  — -  -  — -  (49) 

Pc8  Pc8 

Table  29  illustrates  a  portion  of  the  thermodynamic  properties  of 
the  propellant  system  with  BeH^-  Solution  was  obtained  with  the 

presence  of  three  condensed  phases. 
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